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O ■ ABSTRACT 
(N 

Context. Stellar ages are extremely difficult to determine and often subject to large uncertainties, especially for field low-mass stars. 
■ We plan to carry out a calibration of the decrease in high-energy emissions of low-mass GKM stars with time, and therefore precise 

age determination is a key ingredient. The overall goal of our research is to study the time evolution of these high-energy emissions 
as an essential input to studing exoplanetary atmospheres. 

Aims. We propose to detemine stellar ages with a methodology based on wide binaries. We are interested in systems composed of a 
low-mass star and a white dwarf (WD), where the latter serves as a stellar chronometer for the system. We aim at obtaining reliable 
ages for a sample of late-type stars older than 1 Gyr. 

Methods. We selected a sample of wide binaries composed by a DA type WD and a GKM companion. High signal-to-noise, low- 
resolution spectroscopic observations were obtained for most of the WD members of the sample. Atmospheric parameters were 
determined by fitting the spectroscopic data to appropiate WD models. The total ages of the systems were derived by using cooling 
(~| ' sequences, an initial-final mass relationship and evolutionary tracks, to account for the progenitor life. 

Pi^l Results. The spectroscopic observations have allowed us to determine ages for the binary systems using WDs as cosmochronometers. 

I , We obtained reliable ages for 27 stars between 1 and 5 Gyr, which is a range where age determination becomes difficult for field 

O . objects. Roughly half of these systems have cooling ages that contribute at least 30% the total age. We select those for further study 

since their age estimate should be less prone to systematic errors coming from the initial-final mass relationship. 
^ ' Conclusions. We have determined robust ages for a sizeable sample of GKM stars that can be subsequently used to study the time 
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evolution of their emissions associated to stellar magnetic activity. 

Key words. stars:low-mass, stars: white dwarfs, stars: visual binaries, stars: activity, stars: evolution 
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QQ 1. Introduction modulations). Furthermore, the calibrations are not applicable 

fN| to stars later than K2. 
The age of a star is one of the most difficult stellar parameters to 

determine. Although several age-dating methods ai-e becoming The determination of the a^es of low-mass stars has many 

Q more precise, stellar chronology still holds many uncertainties, applications in astrophysics (Mamajek & Hillenbrand !2008t 

especially when the objects are in the field and moderately old. Soderblom 2010), including the study of galactic populations or 

^ For ages above 0.7 Gyr, open clusters and moving groups are, in ^se of intermediate-age G-, K-, and M- type stars as tracers 

L- general, not very useful as age caUbrators, since they are typi- °^ evolution, such as calibrating of the decrease in high-energy 

cally at large distances, making their individual component stars emissions. Our interest lies in this use of age calibrations for 

difficult to study in detail, paiticulariy ffieir activity characteris- low-mass stars, because of ffie importance of stellar- emissions, 

tics. This is c leariy illustrated in recently proposed calibrations ^nd their time variability, to the understanding of exoplanet at- 

■ (lMamaiek& Hillenbrand 2008) ffiat, while providing good per- mospheres. The host star to a planet is, by far, its main source of 

formance for young ages (15-20% uncertainty), do not yield ac- energy, and many studies (e.g. | Lammer et al. || 2003t | Penz etal. | 

curate results beyond 0.5 Gyr because of the lack of calibrators 2008) have demonstrated that stellar high-energy emissions have 

and the increasing unreliabihty of rotation period estimates (dif- ^ strong impact on planetary atmospheres. Thus, characterizing 

ferential rotation effects or a weak signal from low-amplitude such emissions is centi-al to proper modelling of exoplanet prop- 
erties. High-energy emissions are related to the magnetic activity 
of the star, and this activity decays throughout the stellar lifetime 

~ in an as yet poorly understood manner, making it difficult to cal- 

Based on observations collected at the Centro Astronomico i^jxr^r^um/ •• 
... ,, - ,^,TTAN ^1 Ai J ■ ■ 1 1 1 .J culate the accumulated effects of the UV and X-ray enussion 
Hispano Aleman (CAHA) at Calar Alto, operated ointly by the Max- , , . rr.> i • • • • i > 
Planck Institut fur Astronomie and the Instituto de Astroflsica de °" *e planetary environments. The decay in activity with stel- 
Andaluci'a (CSIC) ^'^^ a ge is intimately lin ked to the rotational evolution of the stars 
** Based on observations made with the WHT (WiUiam Herschel (e.g. Skumanichll972t|^TCs|i997|). Although the time evolution 
Telescope) operated on the island of La Palma by the Isaac Newton of high-energy emissions of solar-type stars is already well con- 
Group in the Spanish Observatorio del Roque de los Muchachos of the strained from The Sun in Time project (Ribas et al. 2005), that of 
Instituto de Astroflsica de Canarias. late GKM stars still needs much improvement. Our objective is 
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to extend The Sun in Time project to cooler stars. This is justified 
by the many differences in the high-energy emissions between, 
e.g., G-type and M-type stars and by the interest of the latter as 
hosts to exoplanets. 

In this work we discuss the first part of our larger plan to de- 
termine age-activity relationships for low-mass stars in general. 
The paper is organised as follows. In §2 we present the results 
for stars with stellar ages below 0.7 Gyr The methodology of the 
age determination for WDs is explained in §3. In §4 we collect 
our sample and describe its selection criteria. In §5 we present 
the observations and the data reduction. Section 6 is devoted to 
the analysis, atmospheric parameter determination, and the sub- 
sequent age determination method. This is followed by §7 where 
we discuss the low-mass companions in the context of defining 
age sequences, and finally in §8 we elaborate on our conclusions. 

2. Preliminary work 

We started by collecting data on the evolution of X-ray emissions 
for the younger stars (<1 Gyr). Their ages were determined from 
cluster and moving group membership (Montes et al. 2001 ), and 
th e logLx values were obtained from a thorough list provided 
bv lPizzolato et al.l (12003). We complemented the measurements 
with X-ray data estimated directly from ROSAT measurements 
for some field stars with age estimates from different methods. 

The young main sequence phase was studied by means 
of stars belonging to the IC2602, IC2391, Pleiades, a Persei, 
Hyades clusters, and the Ursa Majoris moving group. The early 
G-type star sample, as we explained in section 1, was already 
available from The Sun in Time project. These stars have well- 
known rotation periods, temperatures, and metallicities. In gen- 
eral, we only consider stars in narrow spectral type bins (or effec- 
tive temperature) to avoid contamination from the intrinsic vari- 
ation in magnetic activity with stellar mass. The age range <0.07 
Gyr is covered by stars in the IC2602, IC239 1 , and a Persei clus- 
ters. We should point out that these clusters are younger than 
0.07 Gyr, (0.03, 0.03, and 0.05 Gyr, respectively), but since sat- 
uration is present until 0.07 Gyr or longer in the case of K and M 
stars (Stepieri et al. 2001; Jeffries et al. 2010), considering data 
with different ages does not affect the value of log Lx- 

For ages above 1 Gyr cluster or moving group membership 
is not a useful age determination method. Other age indicators, 
such as the use of rotation period, age-activity relations, aster- 
oseismology, or theoretical isochrones, are more useful in this 
age domain. Some of them have been used to obtain the ages 
of a few GKM stars older than ~ 1 Gyr. Only a handful of field 
stars have reliable ages in this domain. These stars are a Cen B 
and Proxima Cen, with ages determined from the isochrone 
and asteroseismologic age of their close companion a Cen A 
(iPorto de Mello et al.ll2008l) . HR7703 shows space motions that 
are typical of a thick disk star, so we can very roughly assume 
an age of 10 Gyr. The X-ray luminosities of these stars were 
de termined from the ROSAT database following the calibration 
in ISchmitt et al.l (Tl 995.) .For stars older than ~6 Gyr long-term 
changes in high-energy emissions are, in general, difficult to dis- 
tinguish from short-term stellar activity variations. 

We have put together all the compilled data, and the evolu- 
tion of logLx with age for three spectral type intervals (GO-5, 
KO-5, and MO-5) is illustrated in Fig. [T] In the case of G type 
stars, the plotted values are quite reliable as they come from a 
thorough analysis of the The Sun in Time sample. For K- and M- 
type stars, however, the plotted values are very crude results cor- 
responding to the few stars described above that just have rough 
age estimates. It is likely that the uncertainty of each point is at 
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Fig. 1. X-ray luminosity versus age for solar-like stars and for 
field and cluster K- and M-type stars. Only a few K and M stars 
are available beyond 0.7 Gyr. 



least of a few tenths of a dex. The figure shows, as expected, that 
M-type stars stay at saturated activity levels for a longer period 
of time than G-type stars. According to these results, solar-like 
GO-5 stars are at saturated emission levels until ages of ~100 
Myr, and their X-ray luminosity decreases rapidly. K-type stars 
have saturated emission levels for a little longer (~200 Myr) and 
then also decrease rapidly. Finally, M0-M5 stars seem to have 
saturated emission levels up to 0.5 Gyr or more and then de- 
crease in an analogous way to G- and K-type stars. 

The interval corresponding to ages youger than ~0.7 Gyr 
is covered well by the used cluster and moving group stars 
(IC2602, IC2391, Pleiades, a Persei, Hyades clusters, and the 
Ursa Majoris moving group), but it is obvious from our pre- 
liminary analysis that a more complete sample of older stars is 
needed to define reliable age-activity relationships. Furthermore, 
it is interesting to note that an age of 0.7 Gyr is a key point for 
our Sun, since it represents the time at which life is supposed to 
have appeared on the Earth's surface. Thus, data for older stars 
of different types will be very important for modelling planetary 
atmospheres in a regime that can be relevant to potential life on 
their surface, as happened to our Earth. 

3. IVIetliodology 

We have developed an age determination method based on the 
use of wide binaries where one of the components, a WD, is used 
as a chronometer The members of a wide binary are assumed to 
have been born simultaneously and with the same chemical com- 
position. Since they are well separated (100-1000 AU), we can 
assume that no interaction has occurr ed between them in the past 
and they have evolved as single stars (Oswa lt et al.ll988l) . We are 
interested in wide binaries composed by a WD and a star with 
GKM spectral type. The evolution of a WD can be described as 
a cooling proc ess, which is relatively well understood at present 
dSalaris et al.l2 000). The total age of the WD can be expressed as 
the sum of its cooling time plus the pre- WD lifetime of its pro- 
genitor Thus, ages can be obtained from an initial-final mass re- 
lationship and stellar tracks to account for the pre- WD lifetime. 
This procedure is analogous to that described by ICatalan et al.l 
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(I2008bl) . It is sound to assume that the age of the WD is the same 
as the that of the low-mass companion, since both members of a 
wide binary were born simultaneously. 

We selected a sample of 30 wide binaries containing a WD 
and a GKM star (see Table [1]). In the sample we favour WD 
components classified as a DA, i.e., with the unique presence of 
Balmer lines in their spectra. As we demonstrate in section 6, the 
fits to the spectral features of these WDs yield realistic values 
for the atmospheric parameters (effective temperature, Teff and 
surface gravity, logg). We have collected a large amount of ob- 
servational data (photometry and spectroscopy) on the WD com- 
ponents to implement the proposed approach. We aim at deter- 
mining total system ages (i.e., cooling ages plus progenitor life- 
time) with precisions of 10-20%. This is sufficient in our context 
since magnetic activity is an intrinsically variable phenomenon 
and any relationship will have some inherent dispersion no mat- 
ter how accurate the ages and coronal/chromospheric fluxes. 



4. Sample selection 

The sample of wide binaries was compilled from the available 
literature. For this we considered the rec ent revision of the N ew 
Luyten Two Tenflis (NLT T) Catalogue (iLuvten 1 979-19801) by 
IChaname & Gouldl mO^, an NLTT catalogue with Hipparcos 
stars from Gould & Chaname (2004) and later completed with 
some pairs from a selection of WDh-M systems made by 
ISilvestiiet ai](l2005h . 

Our first selection was do ne only considering the catalo gues 
of 'Chan ame & Gou ld' (2004) and lGould & Chan ame (2004) and 
following a careful evaluation procedure. First, the WD com- 
ponent had to be classified as a DA (i.e., with the sole pres- 
ence of Balmer lines), so that the fitting procedure can be suf- 
ficiently accurate to derive atmospheric parameters. Secondly, 
the other component of the pair should be a star of spectral 
typ e G, K, or M. Fo r the WD typ e classification we have used 
the iMcCook & SionI (1 19991 l2006h catalogues, and flie spectral 
classification of the low mass stars was made by calculating 
their V-J index and taking the spectral classification into account 
for main sequence stars in Johnson et al.l (ll966l) . Information 
about the two proper motion c omponents of both membe rs of 
the wide binary is availa ble in IChanam6 & Gouldl (I2004h and 
iGould & Chanamel (12004') catalogues for most of the members 
of the sample. We have completed and checked this information 
with data from USNO catalogue and calculated the total proper 
motion of each star. We only consider in our sample the wide 
binaries that satisfy the NLTT proper motion threshold, ju >180 
mas yr ' and that they show the same direction in their proper 
motions. 

The resulting sample was completed considerin g some wide 
binaries from the catalogue of 'S ilvestri et al.l (l2005l) . They stud- 
ied the relationship between age and chromospheric activity for 
189 binary systems composed of a WD and an M-type star For 
this purpose, the authors determined the system's age using WDs 
as stellar clocks b ut assuming a roug h value for their mass, the 
typical O.6M0 ( Silvest ri et al.l I2OOII) . Conversely, the ages es- 
timated are potentially unrehable since the mass of each WD 
is an important parameter for estimating their progenitor life- 
time. With the aim of im proving age est i mates , we performed 
a revision of the sample of ISilvestri et al.l (l2005l) and completed 
our sample by considering the pairs containing a DA WD and 
a GKM star in their catalogue. The targets from Silve stri et alJ 
(I2OO5I) were also screened following the proper motion criteria 
explained in the previous paragraph. 



Table 1. The wide binary sample. Spectral types are estimated 
from the V-J index, except NLTT55287, for which B-V was 
used. 



White dwarf 


V 


Companion 


V 


1/ / 

V-J 




NLTT 1762 


16.59 


NLTT1759 


10.28 


1.46 


Kl 


NLTT10976 


17.20 


NLTT 10977 


13.66 


3.52 


M2 


NLTT13110 


17.23 


NLTT 13 109 


12.63 


3.21 


M2 


NLTT 193 11 


16.58 


NLTT 193 14 


13.18 


3.32 


M2 


NLTT21891 


14.79 


NLTT21892 


15.41 


4.22 


M5 


NLTT26379 


12.92 


NLTT26385 


12.55 


3.75 


M4 


NLTT28470 


13.60 


NLTT28469 


14.09 


4.13 


M5 


NLTT28712 


15.55 


NLTT28711 


15.55 


4.19 


M5 


NLTT28772 


16.67 


NLTT28771 


16.72 


4.18 


M5 


NLTT29967 


17.26 


NLTT29948 


9.96 


1.89 


K4 


NLTT31644 


15.50 


NLTT31647 


16.89 


5.06 


M6 


NLTT31890 


15.86 


NLTT31888 


13.48 


2.95 


MO 


NLTT39605 


16.24 


NLTT39608 


15.74 


3.85 


M4 


NLTT44348 


17.50 


NLTT44344 


11.46 


2.34 


K7 


NLTT56546 


15.90 


NLTT56548 


16.60 


3.99 


M4 


NLTT58107 


16.13 


NLTT58108 


16.13 


5.80: 


M8 


G86-B1B 


16.10 


G86-B1A 


14.26 


3.67 


M4 


LP347-4 


12.92 


LP347-5 


11.70 


3.82 


M3 


LP856-53 


15.00 


LP856-54 


12.18 


2.97 


MO 


LP888-64 


13.56 


LP888-63 


13.80 


4.23 


M5 


WOLF672A 


14.34 


WOLF672B 


14.05 


4.19 


M4 


NLTT4615 


17.48 


NLTT4616 


12.44 


3.13 


Ml 


NLTT7890 


17.39 


NLTT7887 


9.84 


1.90 


K3 


NLTT 15796 


17.21 


NLTT 15797 


15.46 


4.31 


M5 


G107-70 


14.62 


G107-69 


13.52 


4.38 


M6 


NLTT 1374 


16.22 


NLTT 1370 


12.90 


2.19 


K6 


NLTT7051 


16.18 


NLTT7055 


13.06 


2.97 


MO 


NLTT13599 


15.94 


NLTT 13601 


8.42 


2.47 


K7 


NLT55288 


16.50 


NLTT55287 


8.03 


1.21 


G7 


L577-71 


12.80 


L577-72 


13.56 


4.08 


M4 



Our final sample is composed of 30 wide binaries. The sam- 
ple with all relevant information, including the V-J indices of 
the companions and their inferred spectral types, is provided in 
Table [T] For NLTT55287 J photometry was not available and 
the spectral type was determined from the B-V index. In Table[T] 
we provide an approximate value for V-J corresponding to an 
M8 star Although photometric information about the two mem- 
bers of the binary was available in the catalogues, we checked all 
V magnitudes in the SIMBAD database. We found some differ- 
ences or confusions between the magnitudes of the two members 
in some cases. For these targets we made a thorough analysis to 
make sure that each member was correctly identified. Figure 
shows a histogram with the spectral type distribution of the low- 
mass stars. As can be seen, most of the sample is composed by 
M-type stars plus a few K-type stars and one G-type star As we 
discuss in section 6, the WD sample mainly contains DA-type 
WDs, but also a few DC-type WDs. 

5. Observations and data reduction 

We performed optical long-slit low-resolution spectroscopic ob- 
servations for our analysis, with the goal of obtaining measure- 
ments with high signal-to-noise ratio (~75-100). We were able 
to obtain spectroscopic data for 25 WDs in our sample dur- 
ing different observing campaigns between February of 2009 
and February of 2010. The observing runs were carried out 
at Observatorio del Roque de los Muchachos (ORM) and the 
German-Spanish Astronomical Center at Calar Alto (CAHA). 
The 3.5m telescope was used during three nights in February of 
2009 at CAHA, and the William Herschel Telescope (WHT) dur- 
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Table 3. Observation details for the targets in the sample. 



KO - K5 K6 - K8 



X 



Spectral type bins 

Fig. 2. Histogram of the spectral type distribution of the low- 
mass stars in our sample, which are mostly of spectral type M. 

Table 2. Information about the setup used for each observatory. 

Observ. Telesc. Spectrogr. R Arm Grism Spectral 

coverage(A) 

CAHA l5m TWIN 1250 Red T04 5500-7500 

Blue T08 3500-6500 

ORM WHT ISIS 2600 Red R316R 4570-8430 

Blue R300B 2730-6270 



Observed 


Exp. 


S/N 


Num. of 


target 


time (sec) 




exposures 


NLTT1762 


9000 


165 


5 


NLTT4615 


1800 


20 


1 


NLTT7890 


3600 


50 


2 


NLTT 10976 


7200 


75 


4 


NLTT13110 


3600 


80 


2 


NLTT 15796 


7200 


105 


4 


NLTT 193 11 


12600 


160 


7 


NLTT21891 


5100 


320 


5 


NLTT26379 


1200 


330 


5 


NLTT28470 


1500 


290 


5 


NLTT28712 


6300 


165 


5 


NLTT28772 


13500 


190 


9 


NLTT29967 


12600 


185 


7 


NLTT31644 


7200 


285 


6 


NLTT31890 


7200 


140 


4 


NLTT39605 


7200 


120 


4 


NLTT44348 


7200 


100 


4 


NLTT56546 


3600 


140 


2 


NLTT58107 


3600 


190 


2 


G86-B1B 


11700 


210 


8 


G107-70 


3500 


120 


3 


LP347-4 


500 


130 


2 


LP856-53 


4500 


210 


5 


LP888-64 


400 


120 


2 


WOLF672A 


4200 


350 


6 



ing 4.5 nights in February and September of 2009 and February 
of 2010 at ORM. Table|2]provides more details about the obser- 
vation setups used. We managed to observe 25 out of 30 initial 
targets, with the remaining five not visible at the time the observ- 
ing runs were allocated. 

The spectroscopic observations covered the main H Balmer 
lines, from Ha to He or H9, whenever possible. The two spectro- 
graphs used included two separate spectroscopic channels (blue 
and red arms) behind the common entrance slit aperture. We 
chose a proper dichroic position in order to observe the Ha line 
in the red arm and all the other visible Balmer lines in the blue 
arm, covering at least a spectral range of ~3500-7500 A. A suit- 
able grism selection was done to place Ha centred and unvi- 
gnetted for the red arm and a maximum number of unvignetted 
Balmer lines for the blue arm. The slit width for each observa- 
tion was chosen according to the seeing (~l"-2"). Spectra of 
high quality are essential to derive the atmospheric parameters 
with accuracy. We performed as many exposures as necessary 
to guarantee a high signal-to-noise ratio (~75-100) for the final 
spectrum of each object (after the corresponding reduction). The 
exposure times and signal-to-noise ratios for all the targets are 
shown in Table [3] For targets with individual observations with 
low signal-to-noise ratio, we co-added individual 1800 s expo- 
sures to minimise the effects of cosmic ray hits on the CCD. 

The WD spectra were reduced using the standard proce- 
dures within the single-slit tasks of IRAF Q. First, the images 
were bias-corrected and flatfield-normalised, and then, the spec- 
tra were extracted and wavelength-calibrated using arc lamp ob- 
servations. Flux calibrations were not performed since they are 
not necessary for our subsequent analysis, and most of the nights 
were not photometric. We applied the heliocentric correction to 



' IRAF is distributed by the National Optical Astronomy observa- 
tory, which are operated by the Association of Universities for Research 
in Astronomy, Inc., under cooperative agreement with the National 
Science Foundation http : //iraf . noao . edu^ 



each spectrum prior to coadding multiple observations obtained 
in different nights. 

6. White dwarf analysis 

After the corresponding reduction, we carried out a first inspec- 
tion of the spectra. Although all the WDs we selected had previ- 
ously been classified as DA, we found that four of them do not 
show the bluer Balmer lines in their spectra, so they are likely 
to be the DC type. For these four targets we checked the part of 
the spectrum observed with the red arm. The Ha line was not 
visible in any of them, and is faint in DC WDs, the resolution of 
the instrument and the low signal-to-noise ratio of the spectrum 
yields in the absence of such line for these four spectra. Table[T] 
lists the observed 21 DA type WDs in the first set of rows and 
remaining four observed DCs in the second set of rows. 

Once the spectra of the DA targets were coiTectly selected, 
we proceeded to the continuum normalization procedure. The 
normalization has to be carried out carefully because the deter- 
mination of atmospheric parameters is very sensitive to varia- 
tions in the continuum. We carefully defined several clean con- 
tinuum windows by considering the central part of each spec- 
trum and avoiding any irregular feature or the wings of the spec- 
tral lines. We used Legendre polynomials of order 15 to 20 for 
the normalization procedure. 

6.1. Atmospheric parameters 

Before calculating the atmospheric parameters of the WDs (Teft 
and log g), we determined the radial velocities of each star using 
the IRAF task fxcor. Each spectrum was cross-coiTelated with 
a reference model from a giid computed by D. Koester The re- 
sulting radial velocity values can be found in Table|4] The radial 
velocities determined were small in a number of cases (rang- 
ing from 10 to 90 km s ') when compared with the resolution 
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element of our observations. For some of the stars, however, 
the measured radial velocities were higher (ranging from 100 to 
200 km s"') thus becoming relevant to the subsequent analysis. 
In our procedure we took all the radial velocities for consistency 
into account. 

We derived the atmospheric parameters of the WDs by per- 
forming a fit of the observed Balmer lines using the models 
of D. Koester (private communication), following a procedure 
based on;^'^ minimization. The Balmer lines in such WD models 
were calculate d with the modified Stark broadening profiles of 
iTremblav & B ergeron (2009), kindly made available by the au- 
thors. The model grid used covers the range Teff =5 000- 100000 
K, in steps of 1000 K at low temperatures (^eff < 14000 K), steps 
of 2000 K (from = 16000 to 20000 K), steps of 5000 K (from 
Teff = 20000 to 40000 K), and steps of 10000 K for the higher 
temperatures (Tgff > 40000 K). The range of log g is between 5.0 
and 9.0 steps of 0.5 dex. The WD models had been previously 
normalised to the continuum and convolved with a Gaussian in- 
strumental profile with the proper FWHM in order to have the 
same resolution as the observed spectra. The fit of the line pro- 
files to synthetic models was done usin g a subrout ine based on 
the SIMPLEX minimization method o nPress et al.l (i 19921) . We 
observed that the fitting procedure applied to all the observed 
lines (in the red and blue arms) suffered from considerable sys- 
tematic effects. After running a number of tests we decided to 
consider only the lines observed with the blue arm for the analy- 
sis to avoid systematic errors because the Ha line in the red arm 
is less adequate for atmospheric parameter determination. 

Both Teff and log g were obtained simultaneously. It is well 
known that the changes in the profiles of the Balmer lines induce 
a certain degree of degeneracy to the procedure. There are often 
two numerically valid solutions corresponding to minima of the 

function around the maximum strength of the Balmer lines 
(Jeff ~ 12000 K), i.e., a "cool" and a "hot" solution. In some 
cases these two minima are well separated and their individual 
determination is straightforward. But for some other objects, the 
parameters corresponding to the two scenarios are so close as to 
become harder to distinguish. We tried to obtain both solutions 
for the WDs by changing the initial Tetf for the analysis. We 
were able to distinguish between the "cool" and "hot" solutions 
by finding the minimum value of x^ for both cases and com- 
paring the results with the photometric temperatures obtained as 
explained in section 6.2. 

Given that our spectra have a very high signal-to-noise ra- 
tio, the main contributor to the uncertainty in the fitting process 
(besides the adequacy of the models themselves) is the normal- 
ization procedure. The normalization can potentially modify the 
broad profiles of the spectral lines and thus introduce system- 
atic errors in the resulting parameters. We carried out a careful 
analysis to quantify such an error source. Using the IRAF con- 
TiiwuM task, we obtained the blaze functions for each individ- 
ual spectrum. To do so, we devised a procedure that is based on 
normalizing each individual spectrum with the blaze functions 
corresponding to all the remaining spectra taken with the same 
instrument. Each of the resulting spectra was then re-normalised 
(to correct for obvious trends) employing a 5-order Legendre 
function for the observations with higher signal-to-noise and or- 
der 3 for the ones with low signal-to-noise ratio. Finally, each 
of the resulting spectra for every target were fitted to models as 
explained above to determine the atmospheric parameters. The 
final uncertainties were estimated following the prescription of 
[Bergeron et al. ( 1992), i.e., by deriving them from the indepen- 
dent fits of the individual exposures for any given star before 
combination, but in our case by considering all the normalised 



Table 4. Atmospheric parameters determined from fits to the 
Balmer lines to the observed WDs of our sample and radial ve- 
locity values obtained for then. 



WD Name 


Teff (K) 


log g (dex) 


Radial 








velocity (km s"') 


NLTT1762 


10360+200 


8.19+0.15 


185+30 


NLTT13110 


6630±220 


7.99+0.19 


210+15 


NLTT19311 


7600+90 


7.88+0.16 


30+15 


NLTT21891 


13520±580 


7.92+0.07 


110+40 


NLTT26379 


15150±590 


7.91+0.08 


110+40 


NLTT28470 


15290+570 


7.99+0.08 


20+5 


NLTT28712 


9910+250 


8.04+0.16 


125+30 


NLTT28772 


9720+190 


8.09+0.11 


85+20 


NLTT31644 


19180+500 


7.89+0.15 


160+40 


NLTT31890 


10790+290 


8.17+0.07 


5+5 


NLTT39605 


9280+230 


8.12+0.11 


5+5 


NLTT56546 


13660+540 


7.97+0.10 


50+10 


NLTT58107 


11130+380 


8.05+0.09 


50+10 


G86-B1B 


9100+220 


8.12+0.11 


160+20 


LP347-4 


12760+230 


7.91+0.05 


90+30 


LP856-53 


10200+260 


8.26+0.15 


165+30 


LP888-64 


9530+310 


7.92+0.12 


50+25 


WOLF672A 


13330+570 


7.85+0.10 


15+5 


NLTT 10976 


7050+190 


7.68+0.25 


115+20 


NLTT29967 


6180+220 


7.26+0.45 


170+40 


NLTT44348 


6510+190 


7.42+0.33 


-40+10 



Notes. The last three rows correspond to objects for which the fit is 
poor. 



individual spectra. We are aware that this procedure may over- 
estimate the error since some of the renormalised spectra were 
clearly not optimal. However, we prefer to make certain that all 
possible systematic eiTors in the fitting procedure are accounted 
for by taking a conservative approach and considering the un- 
certainties derived in this way. The results of the fits and the 
estimated uncertainties are shown in Table |4] Figures |3] and 
show the fits to the Balmer line profiles for the DA WDs in our 
sample. 

6.2. Photometric effective temperatures 

Some stars in our sample were not visible in the epochs of our 
observing runs. For these stars, and for the ones with no Balmer 
lines in their spectrum, we used the available IR photometry 
to obtain their photometric Teff. The V magnitudes were ob- 
tained from the catalogues used to build the sample and from 
the SIMBAD database, JHK photometry was collected from the 
2MASS All Sky Catalog when available. For V magnitude we 
assumed an uncertainty of 0.02 mag as errors are seldom avail- 
able. 

The Tfiff for each ob ject was obtained by adapting the method 
of Masana et alj (|2006) to WDs. This procedure consists in cal- 
cul ating synthetic photome try using the WD atmospheres of 
Hol berg & BergeronI (l2006h . Subsequently, we developed a fit- 
ting algorithm to compare the observational and theorical values 
and minimise the x^ using the Levenberg-Marquardt method, 
and x^ was defined from the differences between the observed 
and synthetic VJHK magnitudes. The unceitainties in the cal- 
culated photometric temperatures were derived from the covari- 
ance matrix. The VJHK magnitudes and photometric T^ff for the 
sample stars are provided in Table |5] 

In the cases where both spectroscopic data and V and JHK 
2MASS photometry were available, the photometric T^ff was 
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Table 5. VJHK magnitudes and photometric T^ff derived for the stars in our sample. 



WD VJHK T,g{K) 



NLTT 10976 


17.20±0.02 


16.801+0.150 


15.205+0.202 


15.783H 


tO.232 


7220+ 


340 


NLTT19311 


16.58+0.02 


16.070+0.103 


15.833+0.194 


15.583H 


t0.220 


7760+ 


320 


NLTT21891 


14.79+0.02 


14.962+0.041 


15.127+0.055 


15.055H 


t0.167 


12700: 


t670 


NLTT26379 


12.92±0.02 


13.405+0.026 


13.445+0.030 


13.544H 


t0.056 


18880+ 


1210 


NLTT28470 


13.60±0.02 


14.013+0.036 


13.984+0.059 


14.026H 


t0.082 


15640+ 


1150 


NLTT28712 


15.55+0.02 


15.591+0.057 


15.483+0.100 


15.734H 


t0.177 


10700: 


t380 


NLTT29967 


17.26+0.02 


16.018+0.087 


15.684+0.114 


15.500H 


tO.238 


5620+ 


160 


NLTT31890 


15.86+0.02 


15.802+0.067 


15.627+0.136 


15.721H 


t0.217 


9980+ 


320 


NLTT39605 


16.24+0.02 


16.052+0.075 


15.722+0.122 


15.889H 


tO.279 


9100+ 


290 


NLTT56546 


15.90+0.02 


16.571+0.137 


15.995+0.197 


15.522H 


tO.264 


13640+ 


1050 


G86-B1B 


16.10+0.02 


16.038+0.078 


15.565+0.114 


14.562H 


t0.109 


8598+ 


240 


LP347-4 


12.92+0.02 


13.171+0.029 


13.195+0.037 


13.179H 


t0.028 


12640: 


tSOO 


LP856-53 


15.00+0.02 


14.907+0.110 


14.797+0.249 


14.754H 


tO.227 


9800+450 


LP888-64 


13.56+0.02 


13.188+0.067 


13.205+0.061 


13.174+0.098 


8940+ 


190 


WOLF672A 


14.26+0.02 


14.603+0.056 


14.534+0.070 


14.562H 


t0.109 


13600: 


t960 


NLTT4615 


17.48+0.02 


16.466+0.098 


16.245+0.193 


15.555H 


t0.181 


5950+ 


190 


NLTT7890 


17.39+0.02 


16.131+0.081 


16.131+0.081 


16.048H 


t0.165 


5600+ 


160 


NLTT15796 


17.23+0.02 


15.712+0.057 


15.461+0.078 


15.511H 


t0.155 


5025+60 


G 107-70 


14.62+0.02 


13.083+0.022 


12.838+0.022 


12.756H 


t0.025 


5200+60 


NLTT 1374 


16.22+0.02 


16.050+0.082 


15.842+0.168 


15.676H 


tO.248 


9290+ 


320 


NLTT7051 


16.18+0.02 


15.568+0.059 


15.392+0.124 


15.353H 


t0.206 


7630+ 


220 


NLTT13599 


15.94+0.02 


14.598+0.038 


14.232+0.058 


14.136H 


t0.069 


5400+80 


NLTT55288 


16.50+0.02 


15.629+0.067 


15.279+0.086 


15.195H 


t0.168 


6510+ 


170 


L577-71 


12.80+0.02 


13.183+0.033 


12.885+0.033 


12.794H 


t0.034 


11310: 


t300 



Notes. The table is divided into three sections: the first one contains the observed DA WDs, the second one lists the DC WDs, and the third one 
includes the WDs that have no spectroscopic observations. 



G86B1H 

f 


LPa56-D3 




NLTT1762 

fr 


NI.TT10976 

t 


V 

T„=9110 ±220K 
log g=8.1Si0.11 


i- 

T.„=10SOO±360K 
log g = B.S6±0.15 


T.„=9530 +310K 
log e = 7.92±0.1S 


T^=10360±200K 
log g=a. 19+0.15 


i 

T^=7050 ±200K 
log g = 7.6B±0.a5 


NLTT 193 11 

T^=7600 ±100K 
log g = 7,88±0,I6 


NLTT2 1 89 1 

f, 

T.„=135SO±580K 
log g = ?,92±0.07 


Nm26379 

T.„=15150±590K 
log g = 7.91±0.0a 


Nl,TTS8470 

T„=15a90±570K 
log g=7,99±0,08 


NLTT28712 

T.„=9910 ±250K 
log g = 8.04±0.16 


KLTT29967 
V 


NLTT31644 


NI-TT-l'lSda 


NLTT565'16 

|; 


KI.TT5ei07 


T^=61B0 ±320K 
log g = 7.26i0.45 


T„,= 191BO±500K 
log g = 7.89±0.15 


T„,=6l3l0 +190K 
log g = 7.42±0.33 


T^,= l 3660 ± 540 K 
log g = 7.97±0.10 


T„,= 11130±3aOK 
log g = 8. 05+0.09 



A?. (A) 

Fig. 3. Fits to the individual Balmer hnes for the WDs observed 
at the WHT of ORM. The line profiles coiTespond to H/3 (bot- 
tom) up to He, H8, or H9 (top) depending on the object. We have 
applied vertical shifts for clarity. The solid lines are the observed 
spectra and the dashed hnes the model that best fits the profile. 



used to evaluate the reliability of the resulting Teff obtained by 
the line profile analysis. In most of the cases, both Tgff were 
in good agreement, as shown in Fig. |5] NLTT26379 is the ob- 
ject showing the largest deviation. For our subsequent analysis. 




T.„=9110 ±320K 
log g=6.1S ±0.11 



T.„= 12760 ± 230 K 
log g=7.91 ±0.05 




T.„=15150 ± 590K 
log a=7.91 ±0.08 




T„,= 10360 ± 200 K 
log g=8.19 +0 15 



T.,= 7600 ± lOOK 
log g=7 8B +0.16 




T.„=15a90 ± 570K 



T.„=9910 ±a50K 
log g=8.04 +0.16 




Fig. 4. Same as Fig. [3]for the WDs observed at the 3.5-m tele- 
scope at CAHA. 



we consider the Teff value obtained from spectroscopy, since it 
should be more reliable. It is worth noting that the spectroscopic 
analysis of WDs with T^ff < 12000K should be done with special 
care, since their atmospheres could be enriched in helium while 
retaining their DA spectral type (Bergeron et al. 1992). In this 
case, photometric Teff provide the needed cross-check to ensure 
that the overall atmospheric parameters are consistent and accu- 
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2.0x10* - 



^ 1.5x10* - 




I 1.0x10*- y-^fri^ 

/ * 

5.0x10' - - 

'JL \ ^ ^ ^ ^ \ ^ ^ ^ ^ \ ^ ^ ^ ^ \ 

5.0x10-* l.OxlO" 1.5x10" 2.0x10" 

Spectroscopic T^^.f (K) 

Fig. 5. Comparison of photometric and spectroscopic T^s de- 
rived in this work for 15 WDs in our sample. The dashed line 
indicates 1:1 correspondence. The;^'^ of the difference between 
the measurements is 0.57. 



rate. It is reassuring to confirm that our photometric and spectro- 
scopic Teff values are in good concordance for the cooler targets 
of our sample. 

Together with Teff , our procedure requires a value for log 
Given the absence of any other possibility for the objects with 
photometric Teff, we decided to adopt a "mean" logg estimated 
from the WDs in the sample with spectroscopic parameters. We 
employed all WDs with reliable logg values (i.e., the first part 
of TablelHi. We obtained log^ - 8.01+0.20, in which the error 
bar is not the standard deviation but the overall range. 

6.3. Masses, cooling times and total ages 

Once we derived the T^g and log g of each WD, its mass (Mwd) 
and cooling time (fcooi) can be obtained by using the appropri- 
ate cooling sequences. We have adopted the cooling tracks of 
[Salaris et al. (2000), which consider a carbon-oxygen (C/O) core 
WD (with a higher abundance of O at the central core) with 
a thick hydrogen envelope ontop of a helium buffer, q{li) - 
Mh/M = lO-"* and ^(He) = Mne/M = 10"^. These improved 
cooling sequences include an accurate treatment of the crystal- 
lization process of the C/O core, including phase separation upon 
crystallization, together with up-to-date input physics suitable 
for computing WD evolution. The resulting values for Mwd and 
fcooi are listed in Table Q 

Systematic errors related to uncertainties on the CO strat- 
ification and env elope composition may affect the cooling se- 
quences dSalarisI 12009). To check the sensitivity of our results 
to the adopted coolin g tracks, we also used the sequences of 
iFontaine et alJ (l200lh . which consider different core and enve- 
lope compositions. The composition of the cooling sequences 
employed can be found in Table |6l The derived masses do 
not change appreciably when adopting alternative cooling se- 
quences, with differences less than ~0.02 M© over the entire 
mass range. The average of the differences in fcooi were found 
to be around 9%. Thus, we adopted a systematic error of 9% in 
fcooi in our subsequent analysis of the cooling ages. 

For three of our targets, included at the bottom of Table 
our spectroscopic fitting procedure yielded low values of \ogg 
(< 7.7) with relatively poor fits. The photometric temperatures 
available for two of these targets (NLTT10976 and NLTT29967) 
show good agreement with the spectroscopic value, thus indi- 



Table 6. Coohng sequences, and layer composition. 



Sequences 


Core 


H envelope 


He buifer 






?(H) 


qim 


Salaris et al. (2000) 


C/O 


10-" 


10-^ 


Fontaine et al. (2001) 


pure C 


10-" 


10-^ 




50/50 C/O 


10-* 


10-2 




50/50 C/O 


10-'" 


10-2 



eating the consistency of the analysis. Thus, although the logg 
values that we obtained seem reliable, we could not infer reli- 
able m asses for them using the cooling sequences of Salaris et al] 
dloool) since they are out of the covered parameter range. The 
values for Mwd and fcooi obtained via extrapolation for these 
three targets can be found in Table Q We obtained Mwd < 0.5 
Mq, which yields progenitor masses lower than 1 Mq. 

From the derived Mwd and fcooi values, the mass of the pro- 
genitor (Mpio g) can be determin ed using an initial-final mass re- 
lationship (Catalan et al. '2008a'). From the calculated Mpiog, we 
used the stellar tracks of Dominguez et al. (1 19991) to derive the 
progenitor's lifetime (fpiog)- The total age of the WDs, and con- 
sequently that of the low-mass companion, follows directly by 
adding fprog and fcooi- The final age value is provided in Table |7] 
In this Table we compiled the final atmospheric parameters, the 
mass of the WD, the mass of the progenitor star, the cooling time 
of the WD, the main-sequence time of the progenitor star, the to- 
tal age, the ratio between the cooling time and the total age, and 
the V-J colour of the low-mass companions. The first 18 WDs 
in the Table have atmospheric parameters determined from the 
B aimer line profile fitting. The three objects in the middle were 
analysed spectroscopically, but the derived log g are outside of 
the valid parameter range of our procedure. The nine WDs at the 
end of the table have temperatures determined by VJHK pho- 
tometry and the mean log g value for our sample estimated from 
the spectroscopic determinations. 

Some of the WDs in our sample have atmospheric param- 
eters and, in some cases, masses, cooling times, and a ges es- 
timate d from previous i nvesti ga tions, such as Bergero n et all 
(2001), Silvestri etal] (1200 Ih . iGianninas et"al.l (12005), and 
Zhao et al. (201 1). We can compare our results to the previous 
determination. Our Tetf and logg values are compared to those in 
the literature and are shown in Fig. |6] In general the agreement 
between the T^ff determinations is good within the uncertainties, 
but there are some discrepant results for some of the targets , such 
as NLTT26379 and G86-B1B. In the case of G86-B1B, our anal- 
ysis shows good accord between spectroscopic and photomet- 
ric determinations and thus, given the available cross-check, we 
prefer the value we determined. In the case of NLTT26379, our 
photometric and spect roscopic TeifS are dis crepant at the 2.7-cr 
level, and the result of Silvestri et al.l (12001) (using photometry) 
lies squarely in between. Given the better a priori reliability of 
the line profile analysis for Teff above 12000 K, we subsequently 
adopt the spectroscopic T^ff value. The comparison of the log g 
values between our values and the literature yields no very dis- 
crepant results. 

6.4. Photometric distances 

We have collected available parallax measurements for the low- 
mass companions and completed the information with photo- 
metric distances. The photometric distances were determined 
using the evo l utiona ry models for low-mass stars described in 
iBaraffe et al.l (Il997h and using the ages determined for each 
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Table 7. Main results for the WDs in our sample. 



WD 




log,? 


MwD 


^cool 




^prog 


Age 


'cool/Age 


Comp 




(K) 




(Me) 


(Gyr) 


(Mo) 


(Gyr) 


(Gyr) 


(%) 


V-J 


NLTT1762 


10360±200 


8.19+0.15 


0.72+0.07 


0.68+0.10 


2.96+0.52 


0.49+!!!' 

-0.10 


i.2+!;-ii 

-0.14 


58 


1.46 


NLTT13110 


6630+220 


7.99+0.19 


0.59+0.08 


1.61+0.35 


1.64+0.87 


2.3+v-t 


3.9+V-* 

"■^ -1.1 


42 


3.21 


NLTT19311 


7600+100 


7.88+0.16 


0.52+0.07 


0.94+0.13 


0.96+0.73 


12+?n 

-10 


13+?. 


7 


3.32 


NLTT21891 


13520±580 


7.92+0.07 


0.57+0.03 


0.23+0.03 


1.45+0.32 


3.4+!-' 

-1.0 


-1.0 


6 


4.22 


NLTT26379 


15150±590 


7.91+0.08 


0.56+0.03 


0.16+0.02 


1.41+0.36 


3.8+11 

-1.3 


3.9+!-^ 

-1.3 


4 


3.75 


NLTT28470 


15290±570 


7.99+0.08 


0.61+0.03 


0.18+0.02 


1.86+0.38 


-0.29 


1.67+!!-^^ 

-0.29 


11 


4.13 


NLTT28712 


9910±250 


8.04+0.16 


0.63+0.07 


0.62+0.09 


2.05+0.74 


1.15+0" 

-0.34 


1.77+!!-^i 

-0.35 


35 


4.19 


NLTT28772 


9720+190 


8.09+0.11 


0.65+0.05 


0.69+0.08 


2.32+0.55 


0.89+026 

-0.19 


-0.20 


44 


4.18 


NLTT31644 


19180+500 


7.89+0.15 


0.57+0.06 


0.07+0.01 


1.44+0.61 


■^'^-1.7 




2 


5.06 


NLTT31890 


10790+290 


8.17+0.07 


0.71+0.03 


0.60+0.06 


2.88+0.29 


0.53+0!!! 

-0.06 


1 1-1+0.09 
' -0.09 


53 


2.95 


NLTT39605 


9280+230 


8.12+0.11 


0.67+0.05 


0.82+0.10 


2.52+0.55 


0.76+«" 

-0.18 


1 ffC+0.19 

XtJO r\ -,,1 

-0.20 


52 


3.85 


NLTT56546 


13660+540 


7.97+0.10 


0.59+0.04 


0.24+0.03 


1.72+0.45 


1.93+0?' 


2.17+!!-?' 

-0.54 


11 


3.99 


NLTT58107 


11130+380 


8.05+0.09 


0.64+0.04 


0.47+0.05 


2.15+0.45 


1.04+0-30 

-0. Is 


j_5j+0.30 

-0. 19 


31 


5.80: 


G86-B1B 


9110+220 


8.12+0.11 


0.67+0.05 


0.86+0.10 


2.51+0.56 


0.77+0- 18 

-0. 19 


1.64+0-2^ 

-O.ii 


53 


3.67 


LP347-4 


12760+230 


7.91+0.05 


0.56+0.02 


0.27+0.03 


1.36+0.26 


4.3!!:? 


4.5!!:? 


6 


3.82 


LP856-53 


10200+260 


8.26+0.15 


0.76+0.07 


0.80+0.14 


3.26+0.53 


37+0 08 

-0.07 


1 17+0.16 
^•^'-0.16 


68 


2.97 


LP888-64 


9530+310 


7.92+0.12 


0.56+0.05 


0.58+0.07 


1.34+0.52 


4 5+" 


5 r" 

-2.1 


11 


4.23 


WOLF672A 


13330+570 


7.85+0.10 


0.53+0.04 


0.22+0.03 


1.00+0.44 


10!^ 




2 


4.19 


NLTT 10976 


7050+200 


7.68+0.25 


0.40+0.15 


0.99+0.12 


- 


- 


- 


- 


3.52 


NLTT29967 


6180+220 


7.26+0.45 


0.14+0.19 


0.55+0.23 


- 


- 


- 


- 


1.89 


NLTT44348 


6510+190 


7.42+0.33 


0.23+0.14 


0.63+0.20 


- 


- 


- 


- 


2.34 


NLTT4615 


5950+190 


8.01+0.20 


0.59+0.09 


2.40+0.68 


1.73+0.94 


1.9+" 

-0.9 


4.3+7-3 

-1.1 


56 


3.13 


NLTT7890 


5600+160 


8.01+0.20 


0.59+0.09 


2.99+0.78 


1.71+0.94 


2.0+" 

-0.9 


5.0+7-5 

-1.2 


60 


1.90 


G107-70 


5210+ 60 


8.01+0.20 


0.59+0.09 


3.90+1.00 


1.69+0.94 


2.0+2-5 

-1.0 


-1.4 


66 


4.38 


NLTT1374 


9290+330 


8.01+0.20 


0.60+0.09 


0.70+0.12 


1.83+0.92 


-^•"-0.6 


■'••'-0.7 


31 


2.19 


NLTT7051 


7630+220 


8.01+0.20 


0.60+0.09 


1.16+0.22 


1.78+0.93 


-^•'-0.7 




41 


2.97 


NLT13599 


5400+ 80 


8.01+0.20 


0.59+0.09 


3.41+0.87 


1.70+0.95 


■^•"-0.9 


5 4+-« 


63 


2.47 


NLTT15796 


5025+60 


8.01+0.20 


0.59+0.09 


4.60+0.88 


1.68+0.95 


2 1+2-6 


"•'-1.3 


69 


4.31 


NLT55288 


6510+170 


8.01+0.20 


0.60+0.09 


1.77+0.43 


1.75+0.94 


^•''-0.8 


3 6+2-' 

■^•"-0.9 


49 


1.21 


L577-71 


8150+110 


8.01+0.20 


0.60+0.09 


0.98+0.17 


1.79+0.92 


1 6+'-' 




37 


4.08 



Notes. The WDs with r^-ooi/Age ratio greater than 30% are highlighted in boldface. 



tometric temperatures in that the former shows relatively short 
cooling ages (usually < 1 Gyr), while these are significantly 
longer for the latter (~l-5 Gyr). This is not a coincidence since 
older WDs are cooler and the Balmer lines in their spectrum are 
weak or nonexistent. 

To further exploit our sample, it is worth evaluating the re- 
liability of our estimated ages. The uncertainties we provide in 
Table |7] come from the random error estimates related to our fit- 
ting procedure, but they do not include systematic contributions 
potentially arising from the relationships and models used. We 
are not especially concerned with the contribution from mod- 
els since the WD cooling process is relatively simple and main- 
sequence lifetime estimates should be quite reliable. However, 
there is some degree of concern as to the possible existence of 
systematic uncertainties related to the initial-final mass relation- 
ship of WDs. Many improvements have been achieved during 
these last years, for instance with the coverage of the low-mass 
domain (Catalan et al. 2008a). However, there is still a relatively 
large cosmic dispersion in the empirical initial-final mass rela- 
tionship that could be related, e.g., to metallicity effects or the 



wide binary system. The results for the sample are collected in 
Tablel For NLTT193 14 and WOLF672B, we could not directly 
perform a photometric distance determination since their esti- 
mated ages are greater than 8 Gyr, which is the age limit of the 
models we used. Thus, we just adopted a value of 8 Gyr, since 
the evolutionary changes of these stars between 8 and 13 Gyr are 
negligible. 

7. Discussion 

The main product of this part of the work is the determination 
of ages for 27 low-mass stars that are members of wide binary 
pairs. For eighteen of these age determinations we employed 
WD atmospheric parameters coming from a detailled spectro- 
scopic fit, while for nine additional systems we used photometric 
Teff values and an estimated log g from the average of the rest. 
The final age is the sum of two contributions, namely the cooling 
time of the WD and the lifetime of the progenitor, mostly duiing 
the main sequence. As can be seen in Table |7] there is a clear 
distinction between the subsamples with spectroscopic and pho- 
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Comparison of effective temperatures 



• Glonnlnos 


NLTT26379 , 4 , ^ _ 

y 


A Bergeron 




+ Silvestri 




















y \ 


1 .... 1 . . 


5.0x10^ 


1.0x10" 1.5x10" 
T^„ this work (K) 


Comporlson of surfoce grovitles 



• Gianninas 

A Bergeron 
a. 4 - Zhao 




8-0 8.2 
log(g) (this work) 



Fig. 6. Comparison of Teff (top) and log;? (bottoni ) deter - 
mined by us with values obtained by Bergeron et al. ("20011), 
[ Silves tri et al. (2001), Gianninas et al. (2005), and .Zhao et all 
(I2OII . The dashed line is the 1:1 correspondence. 



Table 8. Distances of the low-mass members of the wide binary 
sample. 



Star 


rf(pc) 


NLTT1759 


105±15* 


NLTT 10977 


60±5 


NLTT13109 


20±5* 


NLTT 15797 


40±5* 


NLTT19314 


45±20 


NLTT21892 


30±5* 


NLTT26385 


15+5 


NLTT28469 


15±5 


NLTT28711 


55±15* 


NLTT28771 


50±5 


NLTT29948 


45±5* 


NLTT31647 


25+5 


NLTT31888 


80±5 


NLTT39608 


60±5 


NLTT44344 


55±5 


>JLTT56548 


60±5 



Star 


d (pc) 


NLTT58108 


90+ 


10 


G86-B1A 


25: 


t5 


LP347-5 


8± 


5 


LP856-54 


45: 


t5 


LP888-63 


15: 


t5 


WOLF672B 


15: 


t5 


NLTT4616 


35: 


t5 


NLTT7887 


45: 


t5 


G107-69 


10: 


t5 


NLTT1370 


125+10 


NLTT7055 


65: 


t5 


NLTT13601 


30+ 


5* 


NLTT55287 


35+ 


5* 


L577-72 


15: 


t5 



Notes. * nine low-mass stars with trigonometric parallax information. 



poorly-known mass-loss processes during the AGB phase. For 
this reason, we believe that fcooi is more reliable than fprog. 

Given the considerations above, we favour systems in which 
the contribution of WD cooling time to the total age is at its max- 
imum. In addition, we aim at defining a sample of companion 



Table 9. V-J colour and age for stars with fcooi vs. total age ratio 
>30%, separated in spectral type bins and sorted by age within 
each bin. 



V-J 



Age 



Early K 

1.21 
1.90 


1 n+0.14 
5:0:!! 


Late K 




2.19 


T Q + 1.6 

-Q-7 


2.47 


5 4+- '^ 


Early M 

2.95 

2.97 

2.97 

3.13 


1 1 O+0.08 
^■0-0.8 

4 3+2-3 



V-J 


Age 


Ecirly-ffiid A/ 




3.85 


1 co+0.18 
1.38 „,„ 


3.67 


1 64+''-^'' 


3.21 


■2 Q+2 o 
j.y_, 1 


Mid-late M 

4.18 

4.19 

4.08 

4.38 

4.31 


1 so+0.27 
1.JO_(,20 
1 77+0.71 

26+P-^' 

•'■^-1.4 
6 7+2 ** 



Stars with a wide range of spectral types that cover as much of 
the important 1-6 Gyr age interval as possible. From these two 
premises, we found that a constraint of fcooi/Age > 30% gives 
the best compromise between age robustness and sample size, 
so therefore selected systems fulfilling this requirement. These 
objects are highlighted in boldface in TableQ Eighteen systems 
in our sample fulfil this condition. Almost all of the cooler WDs 
analysed in this work belong to this reliable sample, which is ex- 
pected because their cooling times are relatively long. The last 
column of Table |7]includes, for reference, the V-J index of the 
companion to each WD so that the proper sample selection for 
the next step can be made. Following the same prescription as 
was used in the Sun in Time project, one should define the spec- 
tral type (or mass, or V - 7) intervals and select stars that cover 
a wide age range. Given the V-J colours of our sample stars, 
we grouped them into five bins, namely early K, late K, early M, 
early-mid M, and mid-late M. Table|9]shows these bins with the 
V-J and ages of the stars included in them. Each bin has at least 
two stars (to be added to the cluster averages presented in §2), 
covering roughly the interval between 1 and 5 Gyr These are the 
stars for which we collect activity information so that they can 
be used to trace the time-evolution of their high-energy emis- 
sions from the start of the zero-age main-sequence to at least the 
age of our Sun. This is illustrated in Fig. Q which is the same as 
Fig.[T]but with vertical lines showing the position of the ages that 
will be sampled with the stars we have analysed. The new ages 
cover the relevant area that is currently unpopulated and beyond 
what can be studied with open cluster data. 

8. Conclusions 

We studied a sample of 30 wide binary systems composed of 
a WD component and a low-mass companion. We have ob- 
tained high-signal-to-noise low resolution spectroscopy for 25 
of the WDs, which allowed us to carry out a full analysis of 
their spectra and a fit of WD models to derive atmospheric pa- 
rameters, Teff and \ogg. Total ages for the wide binary systems 
were obtained from a procedure based on the cooling sequences 
of Salaris et al. (2000), the initial-final mas s rel ationship of 
[Catalan et al. (2008a), and the stellar tracks of [Domfnguez et^l 
(1999). We estimate that the total ages for 18 of the systems are 
especially reliable as the relative contribution from the cooling 
time of the WD (which is supposed to be relatively unaffected 
by systematic errors) is greater than 30%. The resulting sample 
covers spectral types from early K to mid-late M and, very im- 
portantly, it contains stars that populate the critical age interval 
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Age (Gyr) 



Fig. 7. logLx vs. age diagram for G0-G5, K0-K5, and M0-M5 
spectral type bins. The vertical lines correspond to the ages of 
the WDs with reliable age determinations from our study. The 
G0-G5 spectral type interval is already fully covered from the 
The Sun in Time project. 

between 1 and 5 Gyr for which the determination of ages using 
other methods is difficult. Most of our targets are within a few 
tens of parsecs, so the study of their activity properties is much 
easier than for stars belonging to intermediate- or old-age clus- 
ters. 

For the selected stars we are obtaining high-resolution spec- 
troscopy to measure chromospheric fluxes and X-ray data so that 
we can buUd reliable relationships describing the time evolu- 
tion of emissions related to stellar magnetic activity. As our fi- 
nal goal, our work should provide critical input information to 
any effort of understanding the long-term evolution of the atmo- 
spheres of planets around low-mass stars, including their poten- 
tial habitability. It will also allow better understanding of the 
evolution of high-energy emissions from low-mass stars. The 
activity-age relationship will allow us not only to address the 
time evolution of high energy radiation, but also to use it as 
an age calibrator when the stellar high-energy emissions can be 
measured. This age indicator is very interesting because it will be 
based on the star's own properties and will be useful in a regime 
(> 1 Gyr) where estimating stellar ages is still very difficult. 
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